Inhomogeneous distribution particles in self-gravitational system. 
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The microcanonical partition function for self-gravitational system in three dimensional case has 
been found. Used approach from the field theory of statistical description of the system was tailored 
to gravitational interacting particles with regard for an arbitrary spatially inhomogeneous particle 
distribution. The entropy of self-gravitational system has been found from extreme condition for 
the effective functional. For inhomogeneous distribution particle (formation few cluster of finite 
size) the entropy are bigger as the entropy homogeneous distribution of particles. The increasing of 
entropy of self-gravitational system after formation cluster motive tendency to disperse. 
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The study of self- gravitational system has as fundamental as practical physical interest. The self-gravitational 
system are interesting for testing ideas about the statistical mechanic description of systems governed by long range 
interaction. A self-gravitational system has also more general problem, studied for a long time 0. The standard 
methods of statistical mechanic cannot be carrier to study gravitational system. Due to this fundamental difference, 
the notion of equilibrium is not always well defined and those system exhibit a nontrivial behaviour with gravitational 
collapse. For system with gravitational interaction the thermodynamically ensemble are inequivalent, negative specific 
heat [H in microcanonical ensemble which not exist in canonical description f3|.In the microcanonical ensemble self 
collapse correspond a " gravithermal catastrophe" and in canonical ensemble to an " isothermal collapse" A self 
gravitational system can increase entropy without bound by developing a dense and hot core surrounded by dilute 
halo. Since equilibrium states are only local entropy maximum. However, if introduce a repulsive potential at short 
distance, complete core collapse is prevented and can to proved that a global entropy maximum now exist for all 
accessible values of energy. The effective repulsion can be introduce in many different way but the physical results 
are rather insensitive to the precise form of regularization. Alternative can consider a classical hard-sphere gas by 
introducing excluded volume around each particle Q. For the gas with pure gravitational interaction between its 
particles the virial coefficients for potential of particles interaction 1/r" if n < 3 and thus the partition function 
diverges, the energy of gravitating particle will not be an extensive parameter. The entropy such system tends to 
infinity too when volume F — > oo. In it was shown that a self-gravitating gas collapses. A nature of the collapse 
and its conditions are explained by using simple and clear consideration The phase transition in such systems 
creates the problem of description in the mean-field thermodynamics approach Two type approaches (statistical 
and thermodynamic) have been develop to determination the equilibrium states of self-gravitational system Q , 0] . 
About all this problem and possible solution very good described in review Q . The collapse in such system are begin 
as spatially homogeneous distribution of particles in the all system at once. Formation of the spatially inhomogeneous 
distribution of interaction particles is a typical problem in condensed matter physics and requires non-conventional 
methods of statistical description of the system was tailored to gravitational interacting particles with regard for an 
arbitrary spatially inhomogeneous particle distribution. This method must employs the procedure to find dominant 
contributions to the partition function and to avoid entropy divergences for infinite system volume. Only a few 
model systems with interaction are known for which the partition function can be exactly evaluated, at least within 
thermodynamic limits Q but not for inhomogeneous distribution of particle. As for existence of equilibrium states, a 
few result have been obtained in the framework of " exact" equilibrium statistic mechanics with one considered a finite 
number of particles [lo| . Is no known exact solution in even in three dimensional self- gravitational system. Formation 
of the spatially inhomogeneous distribution of interaction particles requires a nonconventional method, such as use in 
11 1 , [l3l , ) which are based on Hubbard-Stratonovich representation of statistical sum [iJl ■ This method is now 



extended and applied to a gravitational interacting system to find solution for particles distribution without using 
spatial box restrictions. It is important that this solution has no divergences in thermodynamic limits. For this goal 
can use saddle point approximation which take into account the conservation number of particles in limiting space, 
which provided to nonlinear equation. Described condition give the possibility determine microcanonical partition 
function for self- gravitational system. This partition function in the case homogeneous distribution of particle, and 
in case inhomogeneous distribution with formation few equal cluster with same size has not any peculiarity and fully 
determined thermodynamic parameter of system. Has been shown, that the entropy for inhomogeneous distribution of 
particle in self-gravitational system are bigger as entropy homogeneous distribution of particle. This thermodynamic 
force, which appear by formation spatial inhomogeneous distribution of particle can be motivation of motion of cluster. 
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In particular can consider the microcanonical evolution of self-gravitational system. In the microcanonical ensemble, 
the object of fundamental interest is the density of states, which can written in the standard form [15] 

^E,N = J dqdpSiHip, q) - E)6iNip, q) - N) (1) 

where ^e,n denote microcanonical partition function for fixed energy and number of particles in system. Introduce 
the Lagrange multiplier (3 = ^ inverse temperature and 77 = (3fi where /i-is chemical potential and make Laplase 
transformation can present the microcanonical distribution function in the form : 

nE,N = dp j) dT]exp{PE - T]N} J dqdp exp {-/3i7(p, q) + ryiV(p, q)} (2) 



The entropy and temperature for microcanonical ensemble are defined by S* = In Qe,n and inverse temperature can 

. j_ Aitj jji. v^ao uii. t yjL oj oii^L^L v^tiiixitv^ u-iii- u^ii xv^±±iiula ± — /3 OV E 



determined from relation (3 = The pressure of system defined through formula -P = -Iwp and present the 



equation of state for system. 

A system of interacting particles can be treated in the classical manner as Ising model with Hamiltonian [l^ present 
in term occupation number in the form: 



H{n) 



(3) 



where is the additive part of the energy in the state s which is equal in most cases to the kinetic energy J^] , Wss' 
are interaction energy for the particles in the states s and s' . In this model the macroscopic states of the system are 
described by a set of occupation numbers Ug- Index s labels an individual particle state; and can correspond as well 
as a fixed site on the Ising lattice 17|, which explicit form is irrelevant in the continuum approximation. The number 
of particle is fixed which can determine from relation N{n) — X^s'^s- It is clear that calculating of the partition 
function as a rather complicated problem even in the case of the Ising model. The microcanonical partition function 
of a system of interacting particles can given in the form : 



f2 



E,N 



(b dp j) df] exp {PE - 7]N} exp (-PH{n)) 

{»} 



(4) 



or explicit definition form 



VLe.n ^ f dp i dr/e:cp{PE -r]N}Y^ exp < —P 



(5) 



where ^ implies summation over all probable distributions {ug}. In order to perform a formal summation, additional 
{»} 



field variables can be introduced making use of the theory of Gaussian integrals [ij], [iq : 



exp 



(6) 



where Dip 



and uj„} is the inverse of the interaction matrix. The latter satisfies the condition u;„3,a;s"s' 



Sss' with = ±1 depending on the character of interaction energy. If introduce instead of chemical potential other 
important variable - chemical activity ^ = e^^ — e^, and used the obvious relation dr/ — £,^^d^ the microcanonical 
partition function of a system of interacting particles may be rewritten as : 



nE,N = jdpjdi J D^explpE-N\n^- ^Y.i^^s'^s^s')\l[{^eM-P£s + 'Ps)r' 



(7) 
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This microcanonical partition function can be used for calculating all thermodynamically properties of the system with 
fixed total number of particles and energy of system. Only in this presentation can make summation over occupation 
number. After summation over the occupation numbers [l^ . microcanonical partition function finally reduces to: 



oo 

= ^d/3^de J D^eMPSeffiv.O) (8) 



where can introduce the effective entropy 

Seff - ^ E (^.l'V.¥'.') -Sj2^n{l + ^e^-^^^) + iN+l)ln^-(3E (9) 

s,s' s 

where ^ = e^'^ is absolute chemical activity of chemical potential /i. For this calculation was use the Fermi statis- 
tic, because the two classical particle can not occupied one spatial place. After this we can present partition function 
in the usual form ^e,n = expS', where S is entropy of system. For this goal, for determination microcanonical 
partition function allows use the of efficient methods developed in the quantum field theory without imposing addi- 
tional restrictions of integration over field variables or the perturbation theory. The functional (3Sef / (<p. A) depends 
on distribution of the field variables (p and the absolute chemical activity ^. The field variable ip contains the same 
information as original partition function with summation of over occupation numbers, i.e. all information about 
possible states of the systems. The saddle point method can now be further employed to find the asymptotic value of 
the partition function ^Ie^n for N —^ oo; the dominant contribution is given by the states which satisfy the extreme 
condition for the functional. The particles distribution is determined by the saddle point solutions of equations: 

^^^=0 (10) 
0^ dip 

whether this distribution of particles is spatially inhomogeneous or not. The solutions which correspond the finite 
entropy Seff (<p. A) while the volume of the system tends to infinity, mean such solutions could be thermodynamically 
stable. The above set of equations in principle solves the many-particle problem in thermodynamic limit. The spatially 
inhomogeneous solution of this equations corespondent the distribution of interacting particles. Such inhomogeneous 
behavior is associated with the nature and intensity of interaction. In other words, accumulation of particles in a 
finite spatial region (formation of a cluster) reflects the spatial distribution of the fields and the activity. The inverse 
matrix of the interaction, tOss' = (ks ^ fs' |), in continuum limit should be treated in the operator sense [3], 
i.e. 

^-1 = Srr'L^, = -;r-^A (11) 

where m is the mass of particle and A- Laplace operator. With precision to surface term in continuum case the 
effective entropy takes form: 



Seff{ip,0= J dV I g^J^2^ (Vy>)^ + 6j2^n{l- J^e'^-^^'') \+N\n^~(3E (12) 



As shown before [ll|, in all cases classical Boltzmann statistic for high temperature ^ < 1 and can use expansion 
Y^pln(l+Ce'^-^^p''''^ ^e'^-'^"" + .... Integration over the impulse and coordinates should be performed with 

regard for the cell volume (27r/i)'^ in the phase space of individual states [13] ■ After integration over the impulse the 
effective entropy can present in general form 

Seff O^Jdvl^ 8nGm^i3 ^^^^^ ' ^^^^ | + In ^ - (13) 

where A = ^^p^^ ■ The obtained microcanonical partition function can use to calculation all thermodynamically 

relation for self-gravitational system. At most general case it is impossible, but for individual case can obtain the all 
thermodynamic characteristic of system. Next will be examine separate case. 
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First from all we consider the case of system noninteracting particles, when (p = 0. The effective entropy in this 
case can write in the simple form 

3 

Seff{v,0 = -l dV^l^^y +N\n^-(3E (14) 
The extreme condition = = reduce to two equations 

V^('-^)^=^ (15) 



and 



J 



from which at once determination the chemical activity 

and well-know relation between fixed energy and number of particles in system and inverse temperature 
or ^NkT = E. If determined coefficient substitute in to effective entropy can obtain the ordinary entropy 



, N\ f A'KmE\ ^ 'iN 



(19) 



and the microcanonical partition function for fixed total energy and number of not interacting particles, can write in 
the form 

^^.,. = exp|-ln-(^^j +-| (20) 
which reproduce well-know presentation microcanonical partition function for ideal Boltzmann gas: 

In the case homogeneous distribution particles Vtp = with average distance between them I = can introduce 

the average value of potential 



and determine the ordinary entropy in the form 



The microcanonical partition function can obtain in the form: 
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Microcanonical partition function solve problem determination thermodynamically properties of the self-gravitational 
system in the case homogeneous distribution of particle. From entropy can determine the pressure in self gravitational 
system from homogeneous distribution of particle as Ph = ^^^^^ or equation of state in the form 

O TP 

Pf^Vh = (25) 

But in general case the distribution of particle in self- gravitational system are inhomogeneous. Inhomogeneous 

distribution of particle motivate the long-range gravitational interaction. This system are unstable and all system 
divide as finite state into few cluster finite size. The next task are in developing possible method to determine partition 
function taking into account the inhomogeneous distribution with formation cluster of finite size. Let start with the 
situation when homogeneous distribution particle decomposing to n equal cluster with average volume Vc- Inside this 
volume exist nonhomogcneous distribution of particle with nonzero field variable, and outside this volume the field 
variable are zero because the particle are absence. 

In this case the effective entropy can rewrite in the form: 

Seff{'P,0=n [ ^ dV {-^{^^)^ -^AeA -^A{V ~nVe) + N\n^- 0E (26) 

3 

where = 2'KGm'^P and A = (^^^^ ^ as previously. Minimization of effective entropy on field variables lead to 
next saddle point solutions of equation inside the cluster: 

-^A^ + ^Ae'^ = (27) 



and normalization condition yields: 



n I dV^Ae"^ + ^A{V - nVc) = N (28) 



To multiply the first equation on V<^ and used relation Ay> = V(V<^) can obtain the first integral of this equation in 
the form: 

T^(Vv?)^ + C^e*' = (29) 
where A is unknown integral of "motion" which must determine from physical condition. If use the relation /344 = 



■|j4 the saddle point equation = can rewrite in the form 



d(3 

50 

Va 



n 







dV -^Ae'^ - An + -^A{V - nV,) = HE (30) 



If introduce the density function in the form p{r) = ^Ae'^ can rewrite the saddle point equations in the simple form 
normalization condition 

n j dVp{r) + CA{V - nVc) = N (31) 

and equation for conservation energy 



p{r)-/^'} + ^^A{V-nVe) = PE (32) 



The usual entropy can present in the simple form 

S=^N-2j3E + N\n£, (33) 

For presentation the entropy in term know parameter we must determine unknown reverse temperature /? and chemical 
activity ^. The solution of obtained equation completely solve problem statistical description of self-gravitational 
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system, but in general case this solutions are unknown. To make an attempt to solve this problem in general case was 
take place in article ■ Next will be present a certain solution of self-gravitational system in the case formation few 
cluster same size close pacing particle in ones. This condition correspond the final state of self-gravitational system. 
The field variable inside of cluster are constant and can be present as potential between two close pacing particle 
Lp = ifQ = ^ for r — r' = 2R where R is size of particle. For that is possible used the asymptotic value of field 

variable in center of cluster and determine first integral as = ^^e"^". The initial size of cluster can be determine 
from simple reason. In final case can assume that all particle assemble only to n clusters and take into account the 
finite size of particle which occupied volume Vq = can estimate nVc — NVq. The normalization condition in 

this case can present as 

n^Ae'^° + ^A{V - nVc) = N (34) 

and equation of conservation energy as 

^n^Ae"^" - nVcS.Ae'^" + ^£.A{V - nVc) = E (35) 



from which can obtain the chemical activity 



N 

- V- nV.il -e^o) (36) 



and relation [3E = ^N. The usual entropy can rewrite in the other simple form 



3iV 

S = -N + Nln^- f3E = -N + Nln^- — (37) 
Substitution obtained relation to effective entropy yields: 

N 3N 

S = -N + Nln—- — 7T-^ (38) 

A{V - nVc{l - evo)) 2 ^ ' 

This presentation fully solve problem the statistical description of self gravitational system with formation few cluster 
equal size. Using identity iV — NlnN w InNl as result can obtain the ordinary entropy in the form 

{ 3N 
3Ar;,2(i_iv^(l_,.o))j +2^ (39) 
and the partition function niicrocanonical ensemble for self-gravitational system can present as: 

T/^ 1 AnmEe^ ] ' 



(40) 



If are not gravitational interaction between particles, than e*^" — 1 and the partition function reduces to partition 
function of ideal Boltzmann gas of hard sphere. The equation of state in the case inhomogeneous distribution of 
particle (existence few cluster finite size) can present in the form 

_ 2E 

^"^^^"^^ - 3(i_im(i_e^o)) (^^^ 

After this calculation can contend that the entropy of inhomogeneous distribution of particle (existence few cluster 
finite size) is bigger as entropy homogeneous distribution of particle 

Sink ~ Sh = In TTTT (42) 

(l-^(l-e^o)) ^ ' 

if e'^'* > 1 — -^^^(1 — e'^°) that take place for real self gravitational system. The homogeneous distribution of 
particle to meet the requirements of equilibrium state. This is the thermodynamically reason formation inhomogeneous 



7 



distribution of particle. The relation between ordinary entropy self gravitational system produce the next relation 
between equation of state 

p. , T/. , pVh 

^ mh<'mh ^ (AO] 

If assume that the pressure in both states is equal, come to a determination that the volume of inhomogeneous 
distribution of particle of same self-gravitational system is bigger as volume of homogeneous distribution of particle. 
If the domains is unlimited the density of states diverges when the particles are spread to infinity.. Therefore, there is 
no equilibrium state in strict sense. Self-gravitational system have tendency to disperse. This is already the case for 
ordinary gas in infinite volume. The disperse rate is small in general and the system can be found in quasi equilibrium 
state for a relatively long time. 

Indeed, present equilibrium statistical description tell only dilute structure in the self-gravitational system but not 
describe meta stable states and tell nothing about time scales a kinetic theory. The partition function have not 
any peculiarity for different value of gravitational field. The problem of description of the self-gravitational system 
of particles could be solved with current approach where entropy for finite system could be explicitly calculated. 
Spatial non-uniformity of particles as the equilibrium state might alter necessary activation barrier to proceed with 
transformation when the system is being moved into non-equilibrium state. Gravity factor could cither promote or 
retard such transformation depending on the system and conditions concerned. 
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